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1 Introduction 

One of the intriguing phenomena, observed at HERA, is the behaviour of the energy dependence 
of the ratio jotot in deep inelastic scattering (DIS ). It appears ( see Ref. and Fig. 1 ) 
that this ratio as a function of energy is almost constant for different masses of the diffractively 
produced hadrons over a wide range of photon virtualities . At present there is no theoretical 
explanation for this striking experimental observation. The only valid theoretical idea on the 
market is the quasiclassical gluon field approach (see Ref. 0] ) in which the total as well as the 
diffractive cross section do not depend on energy. 

Recently, Yu. Kovchegov and L. McLerran |^ suggested that the constantcy of the ratio a^^ jotot 
is closely related to strong shadowing corrections ( SC ) for diffractive production. They derived 
a formula for this ratio for the case where only a quark - antiquark pair is produced in diffractive 
dissociation. On the other hand, K. Golec - Bierat and M. Wiisthoff |^ suggested a phenomeno- 
logical model which incorporates two main theoretical ideas regarding the transition between 
"hard" and "soft" processes in QCD 0: (i) the appearence of a new scale which depends on 
energy, and is related to the average transverse monentum of a parton in the parton cascade; and 
(ii) the saturation of the parton density at high energies. This model is successful in describing 
all the available experimental data on total and diffractive cross sections, including the energy 
behaviour oia^^jatot @ i- 

The success of the above mentioned papers prompted us to reexamine the energy behaviour of 
jotot in perturbative QCD ( pQCD ) in more detail. Our approach is based on two main 
ideas used to describe the total and diffractive cross section for DIS: 



1. The final state of the diffractive processes in HERA kinamatic region [0], are desribed 
by the diffraction dissociation of a virtual photon ( 7* ) into a quark-antiquark pair and 
quark-antiquark pair plus one extra gluon (see Fig.2 ); 

2. The Mueller-Glauber approach]^ for calculating SC in the total and diffractive cross sections 
for DIS, this was used successfully to describe other indications of strong SC in HERA data 
0, such as the - behaviour of the Fi slope and energy behaviour of the diffractive 



cross section 10 



The paper is organized as follows. In section 2 we give a simple derivation of the Kovchegov 
and McLerran formula based on the s-channel unitarity constraints. We generalize this formula 
for the case of qqG final state ( see Fig.2 ) and discuss the relation between this approach and 
the AGK cutting rules Section 3 is devoted to a numeric calculation of the ratio jo tot 



without any restriction on the value of produced masses. In section 4 we discuss how limitations 
on the mass range change the energy dependence of the ratio. A summary of our results as well 
as a discussion on future HERA experiments are given in section 5. Appendix fives all formulae 
and all details of our calculations. 
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Figure 1: Experimantal data for the ratio (J^^/(Jtot taken from Ref. j^. 
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Figure 2: Dijfractive dissociation of the virtual photon into quark - antiquark pair ( qq ) and 
quark - antiquark pair plus one extra gluon ( qqG) in pQCD. 

2 Shadowing corrections in QCD 
2.1 Notations and definitions 

In this paper we develop further our approach for diffractive production in DIS started in Ref. 0. 
We use the same notations and definitions as in Ref. 0]. In this paper we examine the main 
physical ideas of Ref. |jl2|, i.e. that the correct degrees of freedom at high enegies ( low x ) are 
colour dipoles, rather than quarks and gluons which appear explicitly in the QCD Lagrangian. 
The consequence of this hypothesis is that a QCD interaction at high energies does not change 
the size and energy of a colour dipole. Hence, the majority of our variables and observables are 
related to the distribution and interaction of the colour dipoles in a hadron. 

To facilitate reading the paper we list the notation and definitions which we will use ( see Fig. 3 ): 



1. denotes the virtuality of the photon in DIS, M the produced mass and W the energy of 
the collision in cm. frame; 

2. xp = [Q"^ + M'^)/W'^ is the fraction of energy carried by the Pomeron ( two gluon ladder 
in Fig. 3 ). Bjorken scaling variable is = Q'^/W^; 

3. We use the symbol x for both xp and xp ■ 

4. f3 = Q"^/ [Q"^ + M^) = xp/xp is the fraction of the Pomeron energy carried by the struck 
quark; 



5. k± denotes the transverse momentum of the quark, and r± = r the transverse distance 
between the quark and the antiquark i.e. the size of the colour dipole; 
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Figure 3: Diffractive production of a quark - antiquark pair. 



6. l± is the transverse momentum of the gluon emitted by the quark ( antiquark ); 

7. z is the fraction of the photon momentum in the laboratory frame carried by the quark or 
antiquark; 

8. bt is the impact parameter of the reaction and is the variable conjugated to q±, the momen- 
tum transfer from the incoming proton to the recoiled proton. Note that t — —qj_; 

9. Our amplitude is normalized so that 



da 



7^|/(s=^y^^)|^ 



(2.1) 



with the optical theorem given by 



47r lmf{s,0) . 
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The scattering amplitude in bt space is defined by 



a''{s,bt) 



1 




f{s,t^-ql) . 



(2.3) 



2ti 



5 



11. The s-channel unitarity constraint then has the form 

2Ima'\s,h) = \a'\s,bt)\' + G'''{s,bt) , (2.4) 

where G^"' denotes the contribution of the all inelastic processes. 
Therefore, in the impact parameter representation: 

atot = 2 J dHt Ima^\sM) ; (2.5) 
(Tel = J dHt\a-\s,bt)\' ; (2.6) 
= / dHtG^^is,bt) . (2.7) 

12. xbG{xb,Q^) is the gluon distribution of the nucleon; 

13. (ydipoie{.xp, r^) is the total cross section for dipole - nucleon scattering, and is given by ( see 
Ref. and references therein ) 

o-dipoie(xp,r^) = -—r\xpG{xp,—), (2.8) 

6 r_|_ 

where xpG(xp, 4-) is the number of gluons with = 4- and energy xp in the nucleon. 

14. To characterize the strength of the colour dipole interaction we introduce 

2 A 

_ (^dipole _ TT as 2 ^ ^DGLAP(^ ^ 

where R is the nonperturbative parameter which is related to the correlation radius of the 
gluons in a hadron. We determine its value from high energy phenomenology and HERA 
experimental data ( see section 3 for discussion ). 

To illustrate the physical meaning of Eq. ( p.9|) we rewrite it in the form: 

Hdipoie = o-Q pir±,x) , (2.10) 

where 

^QDGLAP^^^ 4) 

p{r±,x) = (2.11) 

is the density of colour dipoles of size r± and energy x in the transverse plane. In Eq. (|2.10| ), 
(To denotes the cross section for the interaction of a dipole of size r± with a point - like probe 

TT^as 2 

CO = ^ . 

Hence, Hdipoie is a packing factor for dipoles of size in a proton. If Hdipoie is small, we have 
a diluted gas of colour dipoles in a proton, but at low x the gluon density increases [0 and 
i^dipoie — ^ 1- In a such kinematic region our parton cascade becomes a dense system of 
colour dipoles which should be treated non-perturbatively; 
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15. The amplitude for colour dipole scattering on a nucleon is given by 

O'dipoieis,r±;bt) = ——rlxpG{xp,—;bt), (2.12) 
where xpG{xp, 4-; h) is the number of gluons at fixed impact parameter b^. 



16. It can be shown ( see Ref.|Q and references therein ) that we can write 

a:pG^^^^^(xp,4;&0 = xpG^^^^^(xp,— ) , (2.13) 

r1 r_|_ 



if xpG^'^^^^{xp, bt) satisfies the DGLAP evolution equations [p!^ 



In Eq. ( 2.13|) S{bt) is the nucleon profile function which is a pure nonperturbative ingredient 



in our calculations. We assume a Gaussian form for 



Sih) = ^ , (2.14) 
where R has been discussed above; 

17. For the exchange of one ladder ( "hard" Pomeron ) as shown in Fig. 3, a%poie ^e written 
as 

P I DGLAP 

^dipole = (^dipolei^y^-Lih) = l^dipole {^:''^±)^ (2.15) 

18. (2;, r_L; Q^) is the wave function of the quark - antiquark pair with the transverse distance 
between a quark and an antiquark and with a fraction of energy z ( colour dipole of the 

size ). This wave function depends on the polarization of the virtual photon and it has 
been calculated previously in P W\ . 



^r(;2,r±;Q2) = Qz{\-z)K^{ar^_)- (2.16) 
^J*(z,r^;Q2) = iaK^{ar^y-^^- (2.17) 

where = z{l — z)Q'^ + and subscripts T and L denote the transverse and longitudinal 
polarizations of the photon, respectively. 

19. In our calculations we only require the probabilty to find a quark-antiquark pair with the 
size r_[_ inside a virtual photon, namely 

P'\z,rr,Q') = ^E^JEil^Tp + I^lI'} (2.18) 

'^^ f Ai,A2 

= ^E^?{(^' + (l-^)V^iVr^) +4g^.^(l-.)X(ar^)}. 
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2.2 Shadowing corrections for penetration of qq pair through the tar- 
get. 

2.2.1 General approach 

The physics underlying our approach has been formulated and developed in Refs. [|l^] During 
its passage through the target, the distance r±_ between a quark and an antiquark can vary by an 
amount Ar^ oc Rk±/E, where E is the pair energy and R is the size of the target. Since the 
quark's transverse momentum k± oc l/r±, the relation 

k 1 

Ar^ oc R-^ ^ R < (2.19) 

E r± E 

holds if 

s > 2mR , (2.20) 
where s = W"^ = 2mE with m being the mass of the hadron. 
Eq. ( p. 20 ) can be rewritten in terms of Xp, namely, 

From Eq. ( |2.21|) it follows that rj_ is a good degree of freedom Q for high energy scattering. 

We can therefore write the total cross section for the interaction of the virtual photon with the 
target as follows: 

(^toth*+p) = J J d-'^^±P^'i^^r±'^Q'^) (^dipoleixB,r±) (2.22) 
= 2 I d\ Jdz I £ri_ P^'{z,rr,Q^) Ima%,i,{xB,rrM) ■ (2.23) 
The amplitude for diffractive production of a gg - pair is equal to 

a^^(7*+p -> g + g) = j {z,r^-Q^) ai^^i^{xB.r^-M) "^'Kk^.z.r^) , (2.24) 

where ^'^^(/cx, 2;, r_|_) is the wave function of the quark-antiquark pair with fixed momentum 
and fraction of energy z. To calculate the total cross section of the diffractive production we 
should integrate over all k± and z. Using the completeness of the gg wave function one obtains 

a^^{^*+p^q + q) = J d'h J dz J d' P^\z,rr,Q') \a^,i,{xB,rr,bt)\' . (2.25) 
Utilizing the unitarity constraint we obtain a prediction for the ratio 



^ ^ JdHtJdzJ d'r^ P^'iz,rr,Q') \atoiei^B,rr,bt)\' 
atot 2 / d^ht ! dz ! r_L P^* ((2;, rr, Q^) Im a^J:^p^i^{x b , r±; bt) 
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Eq. ( p.26| ) is a general prediction for the ratio 3? |16||3|. It shows that the diffractive dissociation 



and total cross sections are related through unitarity. However, Eq. ( |2.26| ) is too general to be 
used for pratical estimates. 

Assuming that the dipole - proton amplitude is mainly imaginary at high energy the unitarity 
constraint of Eq. (|2.4|) has a general solution 

el I I \ • / -I n(3:,rj^;fti) 



<,oie{^.r^-M) = i[ \ - e ^ j ; (2.27) 

&Sivoie{^.r^-h) = 1 - e-^(-'^^^^*) ; (2.28) 



where Q is arbitrary real function. 

Using Eq. and Eq. {^^), Eq. reduces to 



^ = = —/ n, TTT • (2-29) 

2 



2 



2 / d%t I dz J d^ {z, rr, Q^) i^l - e 

The main goal of our microscopic approach based on QCD is to find ^2 . 

2.2.2 Mueller-Glauber approach 

One way to get a more detailed picture of the interaction, is to consider the dipole - proton 
interaction in the Eikonal model, which is closely related to Mueller - Glauber approach 

The main assumption of this model is to identify the function Q in Eq. (|2.29|) with the exchange 
of the "hard" Pomeron ( gluon ladder ) given by Eq. ( p.l5| ) ( see Fig. 3 ). Since the gluon 



distribution given by the DGLAP evolution equations originates from inelastic processes of gluon 
emission, we assume an oversimplified structure of the final states in the Eikonal model, namely, 
it consists of only a proton and a quark-antiquark pair ( "elastic" scattering ), and an inelastic 
state with a large number of emitted gluons [Nq oc ln(l/a;) . In particular, we neglect the rich 
structure of the diffraction dissociation processes and simplify them to the final state of p + qq. 
For example, we neglect the diffractive production of an excited nucleon in DIS. 

We will discuss the accuracy of our approach in the next subsection, where we expand our model 
to include an excitation of the target. Our accuracy is restricted by the assumption of only 
a quark - antiquark pair and a nucleon in the final state for diffractive processes. The rough 
estimate for the contribution of all excitations of the nucleon is 



(T^^iY +P qq + N* 



.-(p + p^ P + /V.) ^ ^^^^3 (2.30) 



cr^^(7*+]9— s> qq + N) ^ 2 a^i^p + p p + p) 

consequently, we have to consider the nucleon excitations or, at least, to discuss them in the 
HERA kinematic region. 
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Substituting Q = ( see Eq. ( p.l5| ) ) in Eq. ( 2.29| ) we obtain the Kovchegov - McLerran formula 
|T^, then for large we have 



We can see from Eq. and Eq. ( ^31]) that 



dzPr{z,rr,Q') = X 1. (2.31) 



(7^^ oc R'qy <rl>; 

otot cc R^Q'^I <r'i> ■ 
we can evaluate < r']_> from the condition 

^d^,ole{h = 0) = K^Solfi^^ < ^ ± >) = 1 • (2-32) 

If we assume that xG{x, 4-) depends rather smoothly on r'^, and substitute 4- = as an 
estimate, we have 

^ oc xG{x,Q^) . (2.33) 



< r_L > 

Eq. (|2.33|) gives the ratio 3? constant as well as a^^ oc i^^Q^ x xbG{xb,Q'^)- This simple 



estimate, given in Ref. [§, illustrates that the SC lead to a constant ratio 3?, or, vice versa, the 
constant ratio 3? can be a strong argument for substantial SC. To examine this point we calculate 
the ratio 3? using the same assumption that 4- = in the gluon distribution. Using the result 

of the explicit calculation in Ref. |l^ we obtain: 



^ ^ Of^ ^ ^ _ \n{Kdipoie) + C + {I + l^dipole) Ei{Kdipole) + 1-6 

^toi 2[ln(^^^) + C + (1 + ^i!^)Ei{^^) + 1 - e-^^^] ' 

In Fig. 4 the ratio 3?, given by Eq. ( |2.34|) , is plotted as a function of t^dipoie- One can see that at 
large Hdipoie this ratio has a smooth dependence on i^dipoie- However, the values of i^dipoie that we 
are dealing with are not very large ( t^dipoie ~ 1, see Fig. 5, where t^dipoie is calculated using the 
GRV parameterization for the gluon distribution |]l9l ). Fig.4-a shows that for Hdipoie = 0.2 — 2 
we cannot expect that Eq. ( p. 341 ) to yield a more or less constant ratio 3ft. 



However, we would like to draw the reader's attention to the fact that the ratio differs from the 
small Kdipoie < 1 limit where 3ft oc i^dipoie- 

These simple estimates indicate that the SC are essential, but they are still not sufficiently strong 
to use the asymptotic formulae. We are in the transition region from low density QCD, described 
by the DGLAP evolution, to high density QCD where we can use the quasi-classical gluon field 



approximation pO|. In the transition region the Mueller-Glauber approach is a natural way to 
obtain reliable estimates of the SC. However, we need to study corrections to Eq. ( p.29| ) more 
carefully. The most important among them, are the diffractive production of nucleon excitation 
( see Eq. (|2.3CI| )) and of the qgC system, which gives the dominant contribution to the diffractive 
cross section 0. 
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Figure 4: Ratio cr^^/crtot calculated using Eq. iW^W in the Eikonal model for qq diffractive 
production. 



2.2.3 Diffractive production of nucleon excitations 

A. General approach 

We start with a trivial remark, that the nucleon excitations even in DIS are closely related to long 
distance processes and, therefore, to the "soft" interaction which cannot be determined in QCD. 
An alternate way of saying this is, to attribute these diffractive processes to nonperturbative 
QCD, for which at present we only have a phenomenological approach. Theoretically "soft" 
diffraction can be viewed as a typical quantum mechanical process which occurs, since the 
hadron states are not diagonal with respect to the strong interaction scattering matrix. In other 
words diffractive dissociation occurs, because even at high energy hadrons are not the correct 
degrees of freedom for the strong nonperturbative interaction. Unfortunately, we do not know 
the correct degrees of freedom and below we will discuss some models for them. We denote by 
n the correct degree of freedom or the set of quantum numbers which characterizes the wave 
function These function are diagonal with respect to the strong interaction 

= (^„|T|^'„,) = An Sn,n' , (2.35) 

where parentheses denote all needed integrations and T is the scattering matrix. 

Note, that only for amplitudes An do we have the unitarity constraints in the form of Eq. ( |2.4] ), 
namely, 

/mA:'(.,6,) = \A^:is,b,)\' + Gl\sA) (2.36) 
which has solutions of Eq. (|2.27|) and Eq. (|2.28|) for mainly imaginary An at high energies: 
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Ai =^{l - e-^}; (2.37) 

G';^ = I - e-^"(''^*) . (2.38) 

The wave function of a hadron is 

oo 

^hadron = ^ ^'^ . (2.39) 

n=l 

For a dipole - hadron interaction the wave function is equal to '^dipoie{r±) x "^hadron before colhsion. 
After colhsion the scattering matrix T leads to a new wave function,namely 

oo 

final = ^dipoie(r±) X ^ A„ ^„ . (2.40) 



n=l 



From Eq. ( |2.4(j| ) we obtain the elastic amplitude 



"'dipole 



{'^final\'^dipole{r±) X hadron) = £ tt^ A„(s, 6*) , (2.41) 

n=l 

while for the total cross section of the diffractive nucleon excitations we have 

= {"^ finall"^ finalf " (^/W | ^dipoie (r±) X ^ hadron f ] (2.42) 
oo oo 

= E (^l\Ms,bt)\' - {E alA^{sM)V . (2.43) 

n=l 71=1 

Therefore, using Eq. ( p.31| ) instead of Eq. (|2.29| ), we obtain a generalized formula 

3? = = , 2.44 

2 E;r=i«^ idH,j'i^{i - e — ^} 

One can see that this generalized formula has all the attractive features of Eq. (p.29|) , and at 
small X the ratio tends to 1/2, since the normalization constraint 



E «n = 1 (2.45) 

n=l 

We can estimate using the same Eq. ( p.l5|) here 

= ^^Ipoie e ^ , (2.46) 

where 

<^pole = ^^rlxpG^ixp,^) . (2.47) 
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Therefore, the main difficulty with Eq. ( 2.46|) and Eq. ( 2.47 ) is to determine i?„ and xGn{x, \ 



We only have direct experimental information for the proton radius and the gluon distribution 
in the proton. Unfortunately, we cannot evaluate Eq. ( |2.44| ) without developing some model for 
the diffiactive excitations. 

B. Two channel model for diffractive nucleon excitations. 

The main idea ||2^ of this model is to replace the many final states of diffiactively produced 
hadrons by one state ( effective hadron ). In this case the general Eq. (|2.39| ) reduces to the simple 
form 

^hadron = «! ^^1 + ^2 ^2 , (2.48) 

with the condition a\ + a\ = 1 from Eq. (|2.45| ). 



The wave function of the produced effective hadron is equal to 

= -a2^i + ai^'2, (2.49) 

which is orthogonal to '^hadron- 

Eq. ( p. 441 ) can be rewritten in the form 

DD 

3? = = (2.50) 

J 2 / (r^ ,y,bf) fi? (r^ ,y:bf) \ 

S dH, j '-^ ^a\{l - e-^^Y + al{l - e-^V^fj 

2 JdHtJ ^ ( {1 - } + {1 - e t }j 

with AQ = ^^2 — ^1- 

For Qi and AQ we use the parameterization given by Eq. ( |2.46| ) and Eq. ( p.47| ), and the following 
experimental data and phenomenological observations: 



1. Data for single diffiaction in proton - proton collisions lead to a| ~ 0.2 Ref. 



2. Data for J/\E' photoproduction at HERA |2^ ( see Fig. 6a ) show that t-dependance is 
quite different for elastic and inelastic photoproduction. The measured values are B^i = 
4 GeV-^ while = 1.66 GeV'^. Therefore, we take = 2 B^i = 8 GeV-^ in and 

= 2Bin = 3.32 GeV~^, where R^ is the radius in the exponential parameterization for 
Afi; 

3. The energy behaviour of diffiactive J/\E' photoproduction shows that we can consider this 
process as a typical "hard" process which occurs at short distances [|l^ ; 
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action ( Fig.6-h ). 



4. We use experimental evidence |23| that the cross section for elastic and inelastic diffractive 
J/^' photoproduction are equal. From this fact we can conclude that 



a2(l-«2)(Al])(6t = 0)' i?^ = {l~a2Yni{ht = Q) Ri 



which gives 



m{ht = 0) 
Substituting in Eq. ( p.52|) we find Andipoi 



\ 



■l-ai) fl, 



an 



R 



0). 



D 



3-4 l^dipole- 



(2.51) 



(2.52) 



In Fig. 7 we display the ratio i? as a function of K^ipoie- Comparing this figure with Fig. 4a 
we can conclude that in the two channel model, R increases more or less at the same rate as in 
the Eikonal model, but the value of the ratio depends crucially on the model for the diffractive 
excitation of the nucleon. Fig. 7b shows the contamination of the total diffractive cross section 
by the nucleon excitations. One can conclude that the two channel model gives a small fraction 
of the excitation cross section at sufficiently large Kdipoie- ExperimentallyU],o"f^^j^^^jQ„/(j^^^jj, = 
35 ± 15%. Fig. 7b supports a low value for this ratio, but this question should be reconsidered 
after taking into account the diffractive production of the qqG system. 



II We thank Henry Kowalski for discussing this data and many problems stimulated by these data with us 
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Fig. 7-a Fig. 7-b 

Figure 7: Ratio (T^^/(Ttot (Fig. 7a) and ratio (JeJ^itationl (^elastic (Fig- 7b) versus Kdipoie calculated 
in two channel model for nucleon excitations for qq diffractive production. 

C. Diffractive production in the Additive Quark Model 

As we have mentioned, the main problem of deahng with the "soft" high energy interaction, 
is to find the correct degrees of freedom, and to incorporate them in the general formalism for 
high energy scattering. The two channel model gives an estimate of the importance of proton 
excitation processes, but it is too phenomenological to be instructive and not totally reliable. Here, 
we consider the nucleon excitation in the Additive Quark Model [^. In this model the correct 
degrees of freedom at high energies are the constituent quarks. In spite of a certain naivity this 
model has not been abandoned and it is included in the standard Donnachie - Landshoff Pomeron 
approach for "soft" processes at high energies. 

Inherent in this model is the assumption, that the 7* - constituent quark interactions dominate, 
while other interactions e.g. the interaction of 7* with two constituent quarks, are suppressed by 
factor rq/Rj^ , here rq is the size of the constituent quark and Rn is the radius of the proton. It 
is obvious that in the AQM we have the same Eq. ( p.29| ) ( or Eq. ( |2.26|) ) where Q = Qq describes 
the interaction of a colour dipole with the constituent quark. In the AQM the gluon distribution 
of the constituent quark is equal to xGq{x,Q'^) = ^xG]\f{x,Q'^) and, therefore 

QQ{x,rr,bt) = K%p^i^e = xGn{x,—) e ''Q . (2.53) 

J TT Tq r_|_ 

Consequently, we find rg using the same AQM to describe the double parton cross section mea- 
sured by the CDF |^ (see Fig. 7b ). The CDF collaboration has measured the inclusive cross 
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section for the production of two "hard" pairs of jets, with large and almost compensating trans- 
verse momenta in each pair, and with similar values of rapidity. Such processes cannot occur in 
a one parton shower, and only originate from two parton shower interactions as shown in Fig. 7b. 

The double parton cross section can be written in the form |l27| 



ainei{2jets) ainei{2jets) 
ctdp = rn (2.54) 

where factor m is equal to 2 for different pairs of jets, and to 1 for identical pairs. The experimental 
value|2^ of (Xe// = 14.5 ± 1.7 ± 2.3 mh. 

In the AQM (see Fig. 7b ) (Xe// can be easily calculated and it is equal to 

aeff = 9 X 27rr5 , (2.55) 

where factor 9 reflects the quark counting and 27rrQ comes from the integration over bt. Comparing 
Eq. ( p. 551) with the experimental value of cxe// we obtain Tq = 0.66 ± 0.16 GeV~'^. 

Substituting this result in Eq. ( |2.53[ ) we find n^ipoie ~ 5 n^ipoie "with the same x and r_|_ depen- 
dance. In Fig. 8 one can see the prediction for the ratio a^^/atot- This approach gives the ratio 
which depends smoothly on Kdipoie for Hdipoie > 0.75. 

Comparing Fig. 8 and Fig. 4a we can conclude that the diffractive production of the nucleon exci- 
tations gives about 30% of the total diffractive cross section at Hdipoie ~ 0.2 — 0.6 and becomes 
very small at Hdipoie ~ 1.5 — 2. 



2.2.4 AGK cutting rules and diffractive production 

In this section we derive Eq. ( |2.29| ) in the Mueller-Glauber approach exploiting the AGK cutting 
rules |jll|. This derivation is more complicated than the previous derivation which was based 



directly on the s-channel unitarity constraints. As we intend using the AGK cutting rules for 
calculating the diffraction production of the qqG system, we think it instructive to start with a 
simple example. 

The AGK cutting rules provide a prescription of how to calculate the cross sections for the 
processes with different multiplicities of the produced particles, if one knows the structure of 
the Pomeron exchange, and the expression for the total cross section in terms of multi Pomeron 
exchanges ||2^. In the Mueller-Glauber approach we have such an expression for the total 



colour dipole - proton cross section, namely, (see Fig. 9) 



^dipole 



/„ nP{x.r I -fit) 
d\ { 1 - e ^ } , (2.56) 



where fl is given by Eq. ( |2.15| ). We first need to define the Pomeron structure i.e. to specify the 



kind of inelastic processes that are described by Q^. From Eq. ( p.l5| ) we see that describes 
the inelastic processes with large average multiplicity ( < np > ) of produced partons, since 
it is related to the DGLAP parton cascade. For example, at large and low x, < np > = 
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as InQ^ ln(l/x) ^ 1. The AGK cutting rules allow us to calculate the cross sections for the 
processes with average multiplicities 2< np ?> < np > and so on, as well as the process for the 
diffr active dissociation with multiplicity much smaller that < np >. 

Eq. ( p.56D can be rewritten in the form 

(Td^pole = 2 j d\ i—j , (2.57) 
where each term corresponds to the exchange of n Pomerons. The AGK rules are 

^lr>oieik < np >) = I d\Cni-ir-'^^^^^_ {^''Y ; (2.58) 

(CiP \ " 
}; (2.59) 

where CTdipoiei^ < ^p >) and cr'^j^p^i^^DD) are the contributions of n-Pomeron exchange to the 
cross section for the process with average multiplicity k < np >, and to the cross section for the 
diffractive dissociation processes with small multiplicity. 

Summing over n in Eq. (|2.59| ) we obtain 



. oo / qP \ " 



(2.60) 
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Figure 9: Total cross section for the dipole - nucleon interaction in the Mueller- Glauber approach. 
The dashed line shows the diffraction dissociation cut. 



d%, 2 { 1 - e-T } - { 1 



d\ 1 - e~ — 



(2.61) 
(2.62) 



One can see that Eq. ( p. 621) is just the same equation for diffractive production that we have 
obtained from unitarity ( see Eq. ( p.29| ) ). From Eq. ( 2.58|) we can find the cross section for the 
production of A; - parton showers which is equal to 



CTdipoleik <np>) = j 



kl 



(2.63) 



Eq. ( p.63|) will be very useful below, when we discuss the diffraction production of the qqG system. 



2.3 Cross section for qqG diffractive production with SC 

2.3.1 First correction to the Mueller - Glauber formula for the total cross section 

As shown in Fig. 9, the Eikonal approach takes into account only rescatterings of the fastest colour 
dipole. In this subsection we will extend the formalism so as also to include the rescatterings 
with the target of the two fastest color dipoles: the initial colour dipole and the fastest gluon in 
Fig. 9. Our goal is to include all diagrams of the type shown in Fig. 10. Dashed lines in Fig. 10 
indicate the diffraction dissociation cuts. These diagrams include the diffractive production of 
qqG system as well clS cL quark - antiquark pair. 

Since Eq. ( p^.27[ ) and Eq. ( |2.28| ) give the general solution to the unitarity constraint, our problem 
is to find an expression for fl which will be more general that Eq. ( p.l5|) with /^^glap (jgg^ed by 
Eq. (^. 
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Figure 10: Total cross section for dipole - nucleon interaction in the first iteraction of the Mueller- 
Glauber approach. The dashed lines show the diffraction dissociation cut. 



The natural generalization of 

^DGLAP ig 

to substitute in Eq. (|2.9| ) the Mueller-Glauber formula 
for the gluon structure function |^ namelyQ 



1 ^^DGLAP^^^ ^) ^-1 _^ ^G^/G^^^^.^^) ^ (2.64) 



4 dx' /"^ dr'f ^ , ^ nP(x'y^;bt) 

where 



Jx ~x^ Jrl "Pf ^ 



= W ^) e-^ . (2.65) 



Eq. ( |2.64D takes into account the rescattering of the gluon with the target in the Eikonal approach 
( see Fig. 10), however, the question arises of why we need to only include gluon rescattering for 
the qqG system. To understand the physics of Eq. ( p.64|) it is advantageous to consider the 



equation that describes the rescattering of all partons. Kovchegov p9|] using two principle ideas 



suggested by A. Mueller [12| proved that the GLR nonlinear equation |P is able to describe such 



rescatterings. The principles are: 

1. The QCD interaction at high energy does not change the transverse size of interacting colour 
dipoles, and thus they can be considered as the correct degrees of freedom at high energies; 

2. The process of interaction of a dipole with the target has two clear stages: 

(a) The transition of the dipole into two dipoles, the probability for this is given by 

|*(xoi ^ xo2 + xi2)|' = - , (2.66) 



"It should be stressed that xG{x, bt) is introduced in a such way that / cPbtxG{x, bt) — xG{x, ^). 
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where denotes the size of the dipoles, and z the fraction of energy of the initial 
dipole that the final dipole carries; 



(b) The interaction of each dipole with the target has an amplitude a^'(x, ht, y = ln(l/x)). 



The equation is illustrated in Fig. 11 and it has the following analytic form: 



dadipoiei^oi,bt,y) 
dy 



{2 fldipoie (^02, bt 



2CFas 



TT 



^1 (^dipole{^M,y) + 



Xl2,y) - «dipoie(x02, bt 



xi2,2/) a^-po,e(xi2,bt 



^02 ^12 



:Xo2,2/) } 



(2.67) 



The first term on the r.h.s. of the equation gives the contribution of virtual corrections, which 
appear in the equation due to the normalization of the partonic wave function of the fast colour 
dipole ( see Ref. [0 ). The second term describes the decay of the colour dipole of size xqi into 
two dipoles of sizes xo2 and Xi2, and their interactions with the target in the impulse approxima- 
tion ( notice factor 2 in Eq. (|2.67|) ). The third term corresponds to the simultaneous interaction 



of two produced colour dipoles with the target and describes the Glauber-type corrections for 
scattering of these dipoles. 



=1/ \ 




a (x.,) 



a (x„,) 



Figure 11: Pictorial representation of the nonlinear evolution equation that takes into account the 
rescattering of all partons with the target. 

The initial condition for this equation is given by Eq. ( p.27|) with Q = from Eq. ( p.l5| ) at 

X = Xq. 

In DIS the dominant contribution comes from the decay of a small dipole into two large dipoles. 
Therefore, we can reduce the kernel of Eq. ( p.67| ) to |2^ 



d'-2 — f?^ • (2-68) 

I ^02 ) 



P ^02 ^12 



01 
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We make the first iteration of Eq. ( p.67 ), by substituting tlie Mueller - Glauber formula for color 
dipole rescattering ( see Eq. ( p. 271 ) with Q = from Eq. ( p.l5| ) ), and obtain 



p 



aipoiei f^rst Iteration) = {2 (1 - e"" ) - (1 - e"" )2 } (2.69) 

= Cpas { 1 - } . 



Eq. ( |2.69 ) gives the Mueller - Glauber formula for the gluon structure function of Eq. ( p.64| ) fj\ 



Note, that two assumptions have been made in deriving Eq. ( p.69| ): (i) bt ^ Xi2 or xo2 and (ii) 
we neglected the first term in Eq. (|2.67|) . Both approximations hold in the so called double log 



approximation of pQCD p9[] . 

Eq. (|2.64D describes the passage of the qqG system through the target, as it corresponds to the 



interaction of two colour dipoles, which is our qqG system, but not the gluon. However, in the 
parton cascade in the DIS kinematic region, the gluon always corresponds to two colour dipoles 
of the same size. By changing Q^^'^ we take into account the fact that initial quark 

- antiquark pair can fluctuate in the qqG system many times during the passage through the 
target. This is illustrated in Fig. 10. 

Finally, we obtain the following formula for the total cross section: 

(1) 

where 



'^dipole 



d\2{l - e , (2.70) 



n''''{x,rr,h) = ^rixG^^^(x,ri;6,) (2.71) 
with xG^'^'^{x,r'j_; bt) given in Eq. ( p.64| ). 

2.3.2 Cross section for diffractive production 

We would like to obtain the cross section for the diffractive production of both qq and qqG final 
states using the AGK cutting rules. In Fig. 10 one can see which cuts in the total cross section are 
related to the diffracrive processes. They are shown in Fig. 10 by dashed lines. First, we have to 
generalized the AGK cutting rules, since in Eq. (p.58|) and Eq. ( p.59|) we have used the property 



of one Pomeron exchange i.e. one gluon "ladder" exchange , namely 

2 0^ = Gin , (2.72) 

where G^ stands for the inelastic cross section with large multiplicity < np > { see Fig. 9 ). In 
Eq. ( p.70| ) Q'^'^^ itself has a more complicated structure , which can be recovered using the AGK 
rules of Eq. ( P35| ) and Eq. ( P^PQ . Eq. (^) for fi^^^ gives 

.Sg(.,) = / fi _ e-*V, (2.73) 



vr Jx X' Jrl r'l 



tt 



Eq. ( 2.67 ) is written for large number of colours Nc ^ 1. For finite Nc-, 2 Vl in Eq. ( 2.6£ ) should be replaced 



by 2 nP ^^nP 
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with fig from Eq. (!^M )- 

We denote by fi^f the contribution of all processes given by the AGK rules for Q'^^'^. Using this 
notation the simple generalization of the AGK cutting rules gives the contributions of different 



processes to the total cross section of Eq. (|2.70 



cr. 



(1) 

dipole 



dipole 



0", 



(1) 

dipole 



[n; k) 



>;0) 



n=l 



-1) 



n—k 



{n-k)\k\ 



MG 



(2.74) 
(2.75) 
(2.76) 



where we denote by (J^iloiei^'i ^) contribution of k - cut fi^^*^ to the n-th term of Eq. (|2.74| ), 
cr^^po;g(n; 0) is the cross section for the diffraction production of a quark - antiquark pair for the 
n-th term of Eq. ( p.74|) ( see Fig. 10). The structure of the inelastic processes for each term 
c'^loiei^'^ is rather complicated but is well defined by the AGK rules for fi^^*^. However, we 
need only take the cross section of the diffractive process from each fi*^f or, in other words, we 
should replace ficwf '^Nici^t) in Eq- ( P-75|) . Performing the summation over n and k we obtain: 



DD 



d^poieim q-qG) = I e-^"^(--^*) (e<-(^*) - l) , (2.77) 



where fi*^*^ defined in Eq. ( |2.71| ). Eq. ( |2.771 ) gives the contribution to the diffractive dissociation 
cross section of emission of gluons and k is the number of gluons that we summed over. Eq. ( p.76|) 
leads to 



(^dipoleiQQ ^ QQ) 



(2.78) 



In this paper we are only interested in production of one extra gluon which corresponds to 
Eq. ( ^.75[ ) with k = 1. Therefore 



0" 



DD 

dipole 



{qq - qqG) = J d%, e'^^^^^-^-^*) ^^§(6,) , 



(2.79) 



which we will use for our numerical calculation. 

Finally, collecting Eq. (|]7|), Eq. (|2]7|), and Eq. (|]7|) we obtain the generalization of Eq. (|2]2|), 
which takes into account the diffractive production of both a quark - antiquark pair, and a quark 
- antiquark pair plus an extra gluon final states: 

DD 

3? = ; (2.80) 



CTtot 



a 



DD 



J d\ J d'r^P^\z,rr,Q^)x 



{1 



r + e 



-n''^(x,r^-M) ^2 2a5(r±) dx^ 

Sir Jx X' Jr 



J2 



^ dr^ _t:ci 
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atot = 2 J d\ J d\^P^ (^,r^;Q') ( 1 - e 

We note that Eq. (|2.80|) was derived in the double log approximation to the DGLAP evolution 
equations, in which the colour dipoles in the produced qqG system are much larger than the initial 
quark - antiquark dipole. 

The factor ^-^'''^i^^i'^'M) second term in Eq. (|2.80|) , leads to the suppression of diffractive 



production of the the qqG system. Therefore, in the asymptotic limit at large values of t^dipoie-, 
only elastic rescattering of the quark - antiquark pair survives, leading to the value of the ratio 
jotot — ^ 1/2. It turns out that this factor is important for numrerical calculations in the 
HERA kinematic region ( see the next section). It was omitted in Ref. 0] 



3 Numerical calculation for a^^/crtot 

In this section we present the numerical result for the ratio of the total diffractive dissociation 
cross section to the total cross section in DIS. We postpone to the next section the consideration 
of the influence of the experimental mass cutoff on this ratio. 

The main parameters that determine the value of Kdipoie-, are the value of and the value 
of the gluon distribution. We choose = lOGeV'"^ since it is the value which is obtained 



from "soft" high energy phenomenology [26| |30] and is in agreement with HERA data on J/\E' 



photoproduction For Q^) we use the GRV'94 parameterization and the leading order 



solution of the DGLAP evolution equation |19]. We have two reasons for this choice: 



1. Our goal in this paper is to understand the influence of the SC on the energy behaviour of 
the ratio /atot- However, experience shows that we are able to describe almost all HERA 
data by changing the initial conditions of the DGLAP evolution equations. Unfortunately, 
we have no theoretical restrictions on this input for any of the parameterizations on the 
market. On the other hand, we know theoretically that the SC corrections work in 
a such manner that they alter the initial conditions for the DGLAP evolution, making it 
impossible to apply them at fixed = Qo- With SC we have to solve the DGLAP equations 
starting with Q"^ = Ql{x) where Ql is the solution of the equation i^di^ioif^i^yQoi^)) = 1- 
Therefore, we are in controversial situation: we require a CLAP input but, if we take it 
from the so called global fits, there is a danger that we will incorporate all the effects of the 
SC in the initial condition of these parameterizations. Only data taken after the 1995 runs 
are at energies sufficiently high to effect the low x behaviour of the initial inputs. GRV'94 
is based on the experimental data at rather large x, so we hope that SC are minimal in this 
parameterization ; 

2. The GRV parameterization starts with rather low virtualities ( as low as ~ 0.5GeV~^ 
). This is a major weakness of this approach, since one cannot guarantee that only the 
leading twist contribution is dominant in the DGLAP evolution equations. We agree with 
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this criticism, but the low value of Qq leads to the solution of the DGLAP equation which 
is closer to the leading log approximation, in which we can guarantee the accuracy of our 
master equation ( see Eq. ( p.80| ) ). 

Before we present our numerical results we have to make an important comment. It concerns 
the substitution in Eq. ( |2.64| ), where the iterarted gluon density xG^'-^{x,4:/r'^;b) is defined. 



Due to technical problems, two modifications of the formula are made in the actual numerical 
calculations. First of all we assume the exponential b factorization of xG^'^{x, 6): 



xG^^(x,4/r2;6) = J-xG^'^{xAy)e 



where 

xG'^^ixA/r^) = I d%xG^'^{xA/r^;b) (3.1) 



We justify the factorization by the following two arguments. The first one is numerical. We 
actually have checked that factorization holds numerically with satisfactory accuracy. The second 
argument is that having assumed factorization for xG^^^"^^ on the same basis we can assume 
it for xG^^'^ . The second modification, made with Eq. ( |2.64| ) is due to the fact that the GRV 



parameterization we are working with does not satisfy the DGLAP equation in DLA (see Ref. 
T^ l for the discussion of the problem). The authors of [rH| suggested a modification of Eq. ( |3.1| ), 



which we implement in all our numerical calculations. The final formula for xG^^'^ has the form: 

(3.2) 

TC J X JQq X LJ 

In all the figures representing our numerical results the upper solid line corresponds to the full 
answer, the widely spaced dashed line to diffractive production of qq pair, while the narrowly 
spaced dashed line to the diffractive production of qqG . 

In Fig. [l^ we plot the results of our calculations using Eq. ( p.80| ). We have not added any 
contribution from the nucleon excitations, based on our estimates given in Fig. 7b. In Fig.lTB 



we show our predictions for small values of Q^. The first observation is the fact that the ratio 
increases considerably, reaching the value of about 25 %. It should be stressed that even at rather 
small values of Q^, we do not see any sign of saturation of the energy behaviour of this ratio, 
which is continously increasing in the HERA kinematic region. 

It is interesting to consider separately the quark - antiquark diffractive production and the pro- 
duction of the qqG final state. Both cross sections increase as a function of the energy W. At 
small (Fig. |13[) the total diffractive cross section is dominated (about 70 %) by the diffrac- 
tive production of the quark - antiquark pair. The situation changes when higher values of 
are considered. Then the qqG final state contibutes even more than the qq pair. Indeed, at 
high energy W = 200 GeV, qqG is produced in approximately 25% of all diffractive events at 
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= 8GeV'^. For smaller values of this fraction decreases, being only 10% at = IGeV^. 
This is a expected behaviour if the SC play an important role. The results presented show that 
the contribution of the extra gluon emission is crucial for the predictions, and may lead to a 
hundred percent enhancement of the ratio. 

At small Q^, the fraction of qqG diffractive production to the total cross section increases suffi- 
ciently slow due to the suppression factor e-^" (^''■^'^) in Eq. {^I^. For larger Q^, Q^^^ becomes 



small and the fraction of qqG diffractive production increases faster with the energy. 

In Fig. iMl we illustrate the dependence of our calculation on the value of R^. We would like to 



recall that two values of the radius R, which we use, have the following physics behind them: 

1. R^ = 10GeV~^ is related to the Mueller-Glauber approach, which corresponds to the 
Eikonal model for the nucleon target ; 

2. R^ = 5GeV~^ is the average radius for the two channel model which has been discussed 
in subsection 2.2.3(B) ; 



One can see from Fig. that the value of the ratio /cxtot depends on R^. However, the 



energy dependence is still very pronounced. 

Therefore, the general conclusion of this section is that the SC fail to reproduce the constant ratio 
of cr^^ /atot for DIS seen in the HERA kinematic region. 



4 Ratio (T^^/cTtot in the mass windows 

In this section we examine a possibility that the mass interval will induce an energy independent 
ratio /utot- As one can see in Fig.l, the experimental measurements were made within some 
windows in mass. 

The full derivation of the mass dependent formulae is presented in the Appendix. The cross 
sections for the diffractive dissociation production of qq pair and qqG parton system with a 
definite final state mass are given by Eq. (|A.4|) and Eq. ( [A.10|) respectively. Any mass window 
can be selected for the mass integrals. 

If summation over the whole infinite mass intervals is performed, our formulae for the diffrac- 
tive dissotiation as well as for the total cross sections should in principle reproduce Eq. (|2.80| ). 
However, the two sets of formulas are different but consistent with each other in the leading 
log(l/x) approximation where log(l/a;p) ~ log(l/xB). When xp is replaced by x^, Eq. ( |A.4| ) and 
Eq. ( A. 14 ) reproduce analytically (and numerically) the corresponding expressions in Eq. ( 2.8U| ). 



In the numerical computations we use xp instead of xp since this energy variable reflects the real 
kinematics of the diffraction production process. Since xp = xp/P and the typical values of is 
not very small for qq system, we do not expect large corrections due to this substitution. The qqG 
channel is more sensitive to the kinematic restriction (see Eq. ( |A.10| )). The changes, introduced 
in Eq. ( A. 10 ), concern both the energy variables as well as the integration limits. These cannot 
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Figure 12: The ratio jotot versus W for different values ofQ^ ( = 10 GeV^^ ). The upper 
solid line corresponds to the full answer, the widely spaced dashed line to diffractive production of 
qq pair, while the narrowly dashed line to the diffractive production of qqG . 
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Figure 13: The ratio /(Jtot versus W for small values of Q"^ ( R"^ = lOGeV""^ ) The upper 
solid line corresponds to the full answer, the widely spaced dashed line to diffractive production of 
qq pair, while the narrowly spaced dashed line to the diffractive production of qqG . 



be justified in log(l/x) limit which we used in our general formulae, but we have to introduced 
them to make our calculation reasonable for the diffractive production in the mass window. 



Figure |T5| presents our results for the mass bins, for which experimental data exists (Fig. 1). 
The ratio of the diffractive dissociation cross section to the inclusive cross section is plotted 
as a function of the center of mass energy. The qq pair (tranverse plus longitudinal) and qqG 
contributions are shown separately. 

The ratios obtained do not reproduce the experimental data (Fig. 1). The significant energy 
dependence persists due to the growth of both the qq and qqG contributions. In the wide range 
of the energies our results are smaller than the experimental curves. This observation is quite 
consistent, since our model excluded the target excitations estimated in the Chapter 2 by about 
30%. 

As it should be, at small masses the main contributions come from the qq pair with more than 
50% given by the longitudinal part. The qqG production is suppressed at small masses. Its 
contibution grows with the mass and dominates at large masses. 

Summing all results up to the mass M = 15(GeV) we do not reproduce the inclusive mass 
results of the previous section. This means that even for = 8GeV^ we expect contributions 
from higher masses. At = 8(GeV^) these contributions originated from qqG production. At 
= 60(GeV^) both qq and qqG will be significant above M = 15(GeV). It is seen that only 
about 50% of the inclusive DD production is contributed by small masses up to M = 15(GeV). 

We wish to point that we have completely disregarded possible final state corrections. As an 
example of such corrections, the master formula Eq. ( p.80| ) does not take into account the second 
diagram in Fig. 2. This diagram is believed to give a relatively small contribution compared to 
the first and the third diagrams. 
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Figure 14: Ratto a^^/atot 
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Figure 15: The ratio cr^^/(Jtot o,s a function ofW for differen mass bins ( B? — 10GeV~^ ). The 
upper solid line corresponds to the full answer, the widely dashed line to diffractive production of 
a qq pair, while the narrowly dashed line to the diffractive production of qqG . 



30 



It is worthwhile comparing our model with the Golec-Biernat Wusthoff model 0|, which suc- 
cessefully reproduces the experimental data (Fig. 1). In the Golec-Biernat Wusthoff model the 
effective dipole cross section acwi^,''^), describing the interaction of the qq pair with a nucleon 
has the form: 

(TGw{x,r^) = ao[l - exp(-ri/(4i?2(a;)))]; Ro{x) = (x/a;o)"/' (GeV-^); (4.1) 

(Jo = 23.03 (mb); A = 0.288; Xq = 3.04 ■ 10"^ 

In this model, the diffractive dissociation cross section is given by the squaring of ctcw iii Eq. ( |4.1| ): 

(^Ewi^^r^) = alw/{l6TTBD); Br> = eGeV'^. (4.2) 



A comparison between CW model and the present work model is presented in Fig. |T6. 

We found a significant difference between the two models . The advantage of G-W model is that 
this model takes into account in the simplest way a new scale: saturation momentum Ql{x) = 
4:/Rq{x), but in doing so, this model loses its correspondence with the DGLAP evolution equation. 
Our approach has a correct matching with the DGLAP evolution and we expected that we would 
be able to describe experimental data better than the G-W model. It turns out ( see Fig.l6-a and 
Fig. 16-b ) that we do not reproduce the cr^^ /crtot ratio in contrast to the G-W model, mostly 
due to our 'improvement' in the region of small r^. 

The second remark is the substantial difference in the way we describe the qqG state. We failed 
to find a correspondence between our formula for qqG production which follows from the AGK 
cutting rules, and the G-W description of this process. However, our failure to fit the experimental 
data is mostly due to a large difference in the dipole cross section, rather than in the different 
treatment of the qqG diffractive production. 

Figs. 17-1 and 17-2 show our dipole cross sections at different energies. The main difference with 
the Golec-Biernat and Wusthoff model is the energy rise of the cross sections which follows from 
ht dependance included in the Eikonal formulae but neglected in the Golec-Biernat and Wusthoff 
model. Fig. 17-3 - Fig. 17-10 show whast distances are essential in our calculations. One can 
see that the main contributions in diffractive cross sections stem from longer distances than in 
the total cross sections as have been predicted theoretically 0. As has been known for long 
time the typical distances in diffractive cross section is of the order of the saturation scale Qs{x) 
[K'dipoieijt = '^IQsi.x)) = 1 ) in contrast with total cross section were typical distances are much 
shorter , about 1/Q.. Therefore, one of the reason why we failed to describe the ratio of interest 
could be that our model cannot describe the dipole cross section in the vicinity of the saturation 
scale. However, it was demonstrated in AGL papers ( see Ref. P| ) that our Eikonal model gives 
a good approximation to the correct non-linear evolution equations at that particular distances 
which are essential accordingly to Fig. 17. 



5 Summary and discussions 

In this paper we developed an approach to the diffractive dissociation in DIS , based on three 
main ideas: 
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Figure 16: A comparison between G-W model and our approach: a) The ratio of our calculation of 
the dipole cross sections to G-W model versus r^. b) The same as Fig. 16a but for the diffractive 
dissociation cross sections. All the graphs are plotted for x — 10~^. c) The ratio of the total cross 
sections for Q"^ = SGeV'^ as function of energy W. d) The ratio of the diffractive cross sections 
for — 8GeV^ as a function of energy W. e) and f) the same as c) and d) for — 60GeV^. 



32 



tot 

dipol 



(mb) 




x=10 



x=10 



x=10 



0.1 0.2 0.3 0.4 0.5 0.6 

r (fm) 
Fig.17-1 

Total integrand partial contribution 



3 
2.5 
2 


\ 

\ 


Q^=8 (GeV ^) 
x=10 


1.5 
1 


\ 

\ 




0.5 














1.00 



0.1 0.2 0.3 0.4 0.5 0.6 

r(fm) 
Fig.17-3 

partial contribution 



Total integrand 
30 

25 

20 

15 

10 

5 





\ 

\ 


Q^=8 (GeV ^) 


\ 




\ 




\ / 


-5 

^ x=10 



1.00 



0.1 0.2 0.3 0.4 0.5 0.6 

r (fm) 
Fig.17-5 

partial contribution 



Total integrand 



3 
2.5 

2 
1.5 

1 

0.5 




Q^=60 (GeV ^) 
x=10 




— 1.00 



0.1 0.2 0.3 0.4 0.5 0.6 

r(fm) 
Fig. 17-7 



DD 

dipol 



(mb) 




x=10" 



x=10" 



x=10" 



0.1 0.2 0.3 0.4 0.5 0.6 

r (fm) 
Fig. 17-2 

DD integrand partial contribution 

0.04 



0.03 
0.02 
0.01 




/ 



/ 



\ Q^=8 (GeV ^) 



1.00 



/ 

/ x=10 
0.1 0.2 0.3 0.4 0.5 0.6 

r (fm) 
Fig. 17-4 

partial contribution 
1.00 




0.1 0.2 0.3 0.4 0.5 0.6 

r(fm) 
Fig. 17-6 

DD integrand partial contribution 

0.01 



0.008 
0.006 
0.004 
0.002 




Q^=60 (GeV ^) 




1.00 



0.1 0.2 0.3 0.4 0.5 0.6 

r (fm) 
Fig. 17-8 



Figure 17: o'^ipo/c i Fig. 17-1 ) and ajjlp^i^{rt; x) (Fig. 17-2) versus rt at different values ofx. 

All others figures are the integrants and partial contributions to the cross sections versus distances 
rt at dijfeent values of x and . 
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1. The dominant contribution to diffractive dissociation processes in DIS stems from rather 
short distances [0] and, therefore, we can use pQCD to describe them; 



2. The final state in diffractive dissociation can be simphfied in the HERA kinematic region 
by only considering the production of qq and qqG (0 ); 

3. The shadowing corrections which are essential for the description of the diffractive processes, 
can be taken into account, using Mueller-Glauber approach |Q. 

Using pQCD and Mueller-Glauber approach we derived a generalization of Kovchegov-McLerran 
formula for the ratio cr^^ /(Jtot ( see Eq. ( |2.80| ) ) which is applicable to the HERA experimental 



data on diffractive production at HERA. However, we found that Eq. (|2.80|) cannot describe the 
approximate energy independence of the ratio cr^^ /atot , observed experimentally. Our attempts 
to introduce the experimental cuts for diffractive production does not change this pessimistic 
conclusion. 

Therefore, we believe, this paper is a strong argument that the nonperturbative QCD contribution 
is essential for diffractive production and our approach, based on pQCD, should be reconsidered. 
However, we showed that the main source of the observed energy dependence arises from rather 
short distances ( see Figs. 16 and 17 ) where we did not see any nonperturbative correction to 
the total DIS cross section. In principle, it was pointed out in Ref. that the scale anomaly 
of QCD generates a nonperturbative contribution at high energy at sufficiently short distances 
r1_ ~ 1/Mq ~ 0.25GeV~'^ . We intend studying the influence of such nonperturbative corrections 
in further publications. 

We studied in detail the contribution of the excited hadronic states in diffractive production. We 
found that the contribution of nucleon excitations should depend on and x leading to large 
cross section at bigger values of and higher x. 

We hope that our paper will draw the attention of the high energy community to the beautiful 
experimental data on the energy dependance of the ratio cr^^/atot which have still not recieved 
an adequate theoretical explanation, and which can provide a new insight to the importance of 
nonperturbative corrections at sufficiently short distances. 
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A Appendix 



In this Appendix we present a derivation of formulae for the cross sections for the diffractive 
dissociation production of a qq pair and the qqG parton system when a final state mass window 
is selected. In order to preserve the unitarity relation between the DD cross section and the 
total cross section, we modify the latter. All the formulae are written in the leading log(l/x) 
aproximation of pQCD. Below we only present results derived for the transversely polarized 
photon. The results for the longitudinal part can be obtained by similar treatment. 



A.l qq contribution 

The qq DD cross section has the form [^]: 

- 16 7r 



^ n (2^ L ^vr AT^ (M + g ) X 



M2 + Q2 8 (27r)2 ^ 



5{M^ - 



z{l-z) 

where within our model for the dipole interaction the square of the amplitude is written: 



(A.l) 



2tx 2tx 



(r))*^T (r'){l -e 



n^^'-^(x' ,r: b) 



}{l-e 



ik(f—r') 



(A.2) 



The notation has been introduced in section 2 (see also Fig. 3), while the parameters are defined 
as follows. 



X 



Q^z{l-z) 



It is important to note that in Eq. ( |A.2|) we introduce a correct energy argument x' in VL since 
at fixed the energy of the dipole-proton interaction is (3 W"^. The kinematic constraint forced 
by the delta function sets x' = xp. Performing the angle integrations we obtain 



I/. 



<iq\ 



n'^'"-' jx' ,r;b) 

dr r a Ki{ar) Ji{kr) {\ — e 2 } 



(A.3) 



Finally for the cross-section we get the result 
'a.m. / ^TTT / = / dl? 



a, 



DD 



M3 7o 



with 



AT, 



1 - 4A;2/M2 
: 1 - 2k'^/M'^- 



dr r a Ki{ar) Ji{kr){l — e 



a = Qk/M 



(A.4) 
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A. 2 qqG contribution 

Consider the diffractively produced qqG system with z and z' being the fractions of the energy 
carried by quark and gluon respectively. We assume that the transverse gluon momentum I is 
much smaller than the quark transverse momentum k. In the leading log(l/a;) approximation of 
pQCD z' <^ z. The kinematic constraint is dictated by the final state mass: 

W = k''/z{l-z) + f/z'. 

The cross section for the qqG production is 

with the square of the amplitude 

R is defined as follows. 

= 1^ Tdz' r T^m' - kVz{l -z)- I'/z') e^'-^^-^') X 

TT^ Jo Jr TT TT JO (Stt)"' 

^,{R) ^1{E!) (1 - e-^^) (1 - e-^^) ^J^) as{r'^) (A.7) 



We use the small z' approximation of the gluon wave function 

/o / Dm T3n \ 



Performing the angle integration and removing the delta function by doing the z' integration we 
obtain 

2 1 /• o /-^M^-feVCi-^) dH 

X 



2 1 f f 

" 3M^-kyzil-z)J "^'^Jo (27r) 



2 



J —aKi{ar)^as{r^)rJi{kr)e 2 j^^ ^J2{IR){1 - e 2 ) 



(A.8) 



As a result of the delta function integration we also find that P < 71^2^(1 - 2). It should 
be stressed that r is the size of the initial quark -antiquark pair, while R is the size of the 
produced two colour dipoles. In our approximation, both dipoles have the same size with 3> r. 
Introducing dimensionless variable I 

f - {zM^ - e/{l-z))P (A.9) 
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we finally arrive at the expression for the cross section: 



M^z{l-z) 



dk^ I db^ dl^ X 



(A.IO) 



drr^aK^{ar)^Jas{r^)Ji{kr)e ^ -^UlR){l-e 



Contrary to the qq case, in the present expression x', which is {k"^ / z{l — z) + Q"^) / W"^ is not 
equal to xp. The energy variables are different for Qp^ and ^2^, because they describe different 



physics. Indeed, factor e 



is a probability that qq pair does not interact inelastically before 



emission of the extra gluon, while (1 — e~ 2 y stands for diffractive production of qqG system ( 
two colour dipoles of size R ). For the rescattering of the qq colour dipole the energy is s = W"^ xb, 
while the rescattering of the colour dipoles of size R occurs at energy s' — 13s. 



A. 3 Total cross section 

In order to preserve the unitarity relation between qq DD cross section and the total cross section 
we modify the latter. Similarly to the qq case we write for the total cross section 



HM^ - jT^) \W (A.11; 



ik{f—r') 



with the square of the amplitude 

n r n f d'^r d'^r' . . , n-^^Gj^/ ^ n'^^^a,' ,T';b) 

\ftot? = / ^ V ^ ^-^2 " ' ' 

(A.12) 

j d% drrdr'r' Ki{ar) Ji{kr)Ki{ar ) Ji{kr'){2 - e ^ - e ^ } 

One of the space integrations can be performed analytically noting that 

/ drrKi{ar)Ji{kr) = , (A.13) 

^0 a [k^ + a^) 

Substituting a = Q k / M we finally obtain the total cross section 

fdM Q /-^/^s „ 1 - 2k^ / 

atot — 16 Ctem / -TT^ / k dk . X 

emj ^2 Q2 + ^^^2 ^l-Ak^/M^ 

/„ r n'^'^ix p. r-.b) 
db^ / drrKi{ar) Ji{kr) {1 - e (A.14) 
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In the above expression for the total cross section the mass integration should be carried out over 
the whole infinite mass interval. The result obtained is consistent with the previous expression 
for the total cross section. If we change xp to xb and perform the M integration we reproduce 
the old result written in Eq. ( p.SOp . A useful equality we use is 

/ n2 + '^'^^''^ = QK,{Qr) (A.15) 
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